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摘 要 : ÆT Google Earth Engine(GCEE) 遥 感 云 平 台 ,利用 2000 一 2019 年 MODIS 积 雪 产 品 资料 提取 
和 计算 新 疆 积 雪 终 日 信息 ,利用 趋势 分 析 , 变 异 系数 等 方法 分 析 了 新 疆 积 雪 终 日 时 空 变化 特征 和 
变化 趋势 。 结 果 表 明 :(1) 新 疆 积 雪 终日 以 天 山 为 界 , 天 山 以 北 长 于 南部 ,山区 为 积 雪 终 日 的 高 值 


区 ,盆地 为 积 雪 终 日 的 低 值 区 。 北 履 准 噶 尔 盆地 和 伊犁 河谷 积 雪 终 日 在 73~114 d 之 间 , 南 绪 塔 里 
木 盆 地 在 0~31 d 之 间 属 于 低 值 区 。 阿 尔 泰 山脉 .天 山 山 脉 和 昆仑 山脉 区 域 在 224~365 d 之 间 属 于 
高 值 区 。(2) 南 强 和 和 北 疆 积 雪 终 日 有 明显 的 时 空 差异 ,2000 一 2019 年 北 疆 准噶尔 盆地 和 和 高 海拔 山脉 
地 区 积 雪 终 日 有 明显 的 推迟 趋势 ,推迟 幅度 达到 14d, 占 新 疆 总 包 
疆 区 域 有 明显 的 提前 趋势 ,提前 幅度 达到 16d 约 占 新 疆 总 


积 的 8%。 南 疆 塔 里 木 倪 地 和 和 东 
积 的 44%。 塔 里 木 盆地 和 准噶尔 盆地 


具有 相反 的 变化 趋势 (3) 新 疆 积 雪 终 日 年 际 变化 差异 显著 ,天 山中 上 段 和 北 疆 积 雪 终 日 出 现 不 稳定 
状况 ,天 山中 段 2002 一 2009 年 总 体 上 呈现 “5M” 型 的 特点 , 即 多 年 积 雪 消 融 日 年 均值 中 出 现 明显 的 
波峰 和 波 谷 , 北 疆 2009 一 2019 年 积 雪 终日 有 较 大 的 年 际 变化 呈现 出 不 稳定 状况 ,出 现 明显 的 波峰 


和 波 谷 ,年 际 变化 较 大 。 
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文章 编号 : 


积 雪 是 冰冻 圈 的 主要 组 成 部 分 ,具有 较 大 的 年 
际 变化 。 季 节 性 融雪 的 变化 对 世界 的 气候 水 文 、 
生态 和 人 类 进程 产生 了 相当 大 的 影响 。 年 度 季 闻 
性 融雪 的 时 间 以 及 该 时 间 的 任何 潜在 变化 都 具有 
广泛 的 生态 意义 ,因此 会 影响 人 类 的 生计 。 在 高 纬 
度 和 高 山区 ,在 水 文 循环 最 重要 的 阶段 中 ,冬季 积 
雪 的 年 融 水 为 全 球 约 1/6 人口 提供 了 主要 的 水 流 和 
地 下 水 补给 源 。 北 极 和 南极 融雪 的 时 间 影 响 了 
生长 期 ,而 且 整 个 冬季 持续 的 积 雪 使 植被 免 受 恶劣 
温度 和 风 的 影响 “。 这 种 季节 性 的 水 资源 是 在 地 
球 变 暧 的 气候 下 变化 最 快 的 水 文系 统 之 一 ,将 广泛 
影响 区 域 经 济 、 生 态 系 统 功能 并 增加 潜在 的 洪灾 人 危 
害 ”““。 预 期 较 高 气温 将 改变 降水 的 类 型 和 开始 时 
间 ; 包 括 北美 落 基山 脉 在 内 的 多 个 地 区 已 经 测算 出 
积 雪 量 减少 ,而 且 气 温 升 高 使 降水 从 降雪 转向 降 
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雨 ,导致 融雪 更 早 发 生 ”。 冰 雪 融 水 是 新 疆 干 旱地 
区 农业 的 命脉 ,区 域 气候 变化 引起 的 融雪 变化 对 春 
夏季 河川 径流 产生 影响 ,将 对 干旱 地 区 的 经 济 和 及 
弱 的 生态 环境 造成 严重 后 果 , 其 至 导致 旱 小 灾害 频 
发 ""。 春 季 融 雪 时 间 异 常 正在 对 包括 农业 在 内 的 
水 资源 可 持续 性 造成 严重 威胁 ,其 至 会 引起 水 文 
极端 事件 ,比如 春季 洪水 ”。 每 年 季节 性 融雪 的 时 
间 以 及 时 间 上 的 任何 潜在 变化 都 具有 广泛 的 生态 
影响 ,因此 会 影响 人 类 的 生计 ,尤其 是 在 高 纬度 地 
区 和 山区 系统 及 其 周围 地 区 。 高 山 雪 融 水 是 新 疆 
内 陆 河 流水 资源 的 主要 来 源 “。 因 此 , 近 些 年 来 ， 
对 新 疆 积 雪 的 研究 较 多 ,相关 研究 结果 显示 新 疆 地 
区 积 雪 最 大 深度 呈现 增加 趋势 ,而 积 雪 覆盖 期 变 
短 “” ,但 大 都 主要 针对 新 疆 局 部 区 域 ,如 北 疆 、 天 
山 山 区 ”阿勒泰 ”等 进行 探讨 和 研究 。 易 颖 等 ” 
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使 用 MODI 数据 产品 提取 了 叶 尔 芜 河流 域 积 雪 物 
候 参 数 ,系统 地 分 析 了 其 空间 分 布 与 变化 特征 。 梁 
鹏 斌 等 2 利用 遥感 数据 资料 和 和 气象 数据 ,分 析 祁 连 
山 积 雪 面 积 动态 变化 特征 ,并 与 气温 降水 关系 进行 
了 分 析 和 研究 。 目 前 大 多 数 国内 外 人 研究 主要 聚焦 
于 流域 尺度 的 积 雪 特性 ,融雪 产 流 过 程 及 融雪 径流 
模拟 等 ”下 。 针 对 新 疆 积 雪 消 融 过 程 的 研究 ,前 人 
多 采用 站 点 积 雪 数据 与 站 点 气象 数据 (点 -点 ) 的 方 
法 开展 新 疆 融 雪 变 化 及 径流 模拟 研究 ,大 多 数 基于 
流域 尺度 的 以 雪 特 性 .融雪 产 流 过 程 及 融雪 径流 模 
拟 等 ,并 取得 了 一 定 的 成 果 。 然 而 针对 新 疆 积 雪 终 
日 利用 遥感 技术 研究 较 少 ,以 年 为 时 间 尺 度 的 长 时 
间 序 列 的 积 雪 终 日 时 空 变化 趋势 研究 则 更 少 。 研 
究 新 疆 积 雪 终 日 对 气候 变 暖 的 响应 ,不 仅 对 亚洲 区 
域 和 全 球 气 候 变 化 的 诊断 具有 重要 意义 ,也 能 为 流 
域 下 游 的 农业 ,生态 以 及 防 灾 减灾 等 提供 信息 服务 。 

本 文 基于 GEE 的 云 计算 平台 ,利用 MODIS 每 日 
积 雪 产品 ,像素 为 单位 计算 每 年 的 无 雪 日 ,跟踪 积 
雪 终 日 年 际 变化 ,在 长 时 间 尺 度 上 分 析 新 疆 积 雪 终 
日 的 时 空 变化 ,探讨 近 20 a 全 球 气候 变 暧 背景 下 新 
吐 各 地 积 雪 终日 的 时 空 变化 规律 ,从 而 更 好 地 了 解 
新 疆 积 雪 终日 时 空 变化 如 何 响应 气候 变化 ,为 应 对 
气候 变化 以 及 融雪 水 资源 的 优化 管理 提供 理论 
依据 。 
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sj 


1 研究 区 概况 


新 疆 地 处 亚 欧 大 陆 中 心 .中 国 西北 边 隆 ,地 处 
73°45’ ~96°39' E, 34934'~49s17'N ,总 面积 约 1.66x 
10'km。 新 疆 北 为 阿尔 泰山 脉 , 南 为 昆仑 山脉 ,天山 
山脉 横 吾 中 部 ,把 新 疆 分 隔 为 北野 和 南 疆 两 大 区 
域 , 北 狂 的 准噶尔 盆地 和 南 狂 的 塔里木 盆地 与 此 三 
座高 山 构 成 了 “三 山 夹 两 使 ”的 独特 地 貌 格局 (图 
1)。 新 疆 具 有 典型 的 内 陆 干 旱 和 半 王 旱 气候 特征 ， 
冰川 积 雪 才 盖 的 山脉 环绕 在 盆地 的 边缘 ,而 沙漠 和 
湖泊 则 在 盆地 中 广泛 分 布 。 新 疆 隆冬 时 期 积 雪 履 
盖 面 积 为 100x10' km ,冬季 雪 贮 量 在 全 国 最 为 丰富 
达 181.8xl08m ,水 当量 占 全 国 冬季 平均 积 雪 贮 : 量 的 
33.9%'”。 全 疆 降 水 量 的 时 空 分 布 不 均匀 日 年 均 降 
水 量 较 少 , 约 145 mm, 其 中 北 疆 约 150~200 mm, 南 
疆 不 到 100 mm, 26 ŒX 2000~2500 mm ,干旱 度 为 
4~16"”"。 由 于 特殊 的 地 理 位 置 .地 形 和 气候 条 件 ， 
新 疆 总 体 上 生态 系统 比较 脆弱 敏感。 


2 数据 与 方法 


2.1 数据 来 源 
(1) 积 雪 数据 :本 文 使 用 的 MODIS 全 球 每 日 积 
雪 产 品 数据 来 自 于 美国 国家 冰雪 数据 中 心 (https:/ 


图 号 为 新 S(2021)046 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 


图 1 研究 区 示意 图 
Fig. 1 Sketch map of the research area 
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nsidc.org/data/MOD10A1 ) 的 MODIS/Terra 合成 积 雪 
产品 MODIS10A1 ,空间 分 辨 率 为 500 m; 时 间 范 围 为 
20004F 1A 1 A—20194F 12 H31 H. RA iE 
使 用 归 一 化 积 雪 指数 (Normalized difference snow in- 
dex, NDSI) 和 一 系列 旨 在 减轻 错误 并 标记 不 确定 的 
队 雪 检测 结果 来 识别 积 雪 。 积 雪 终 日 为 每 年 最 后 
一 次 出 现 积 雪 的 日 期 , 积 雪 终 日 是 指 陆 面 上 无 任何 
像 元 被 识别 为 雪 履 盖 的 日 期 ,文中 将 当年 1 月 1 日 
开始 至 12 月 31 日 结束 定义 为 一 个 积 雪人 年。 

(2) DEM 数 据 :DEM 数据 由 美国 国家 航空 航天 
局 (NASA) 与 日 本 经 济 产业 省 (METD 合 作 完 成 的 先进 
星 载 热 发 射 和 反射 辐射 仪 全 球 数 字 高 程 模型 AS- 
TER GDEM V2.0( 分 辩 率 为 30 m) 数 据 下 载 获取 。 对 
高 程 数据 进 了 投影 转换 , 重 采 样 等 处 理 ,最 后 统一 
为 500 m 空 间 分 辨 率 ,得 到 与 积 雪 数据 分 辩 率 和 投 
影 系统 相同 的 空间 数据 。 
2.2 研究 方法 

本 文 积 雪 终日 参数 是 使 用 平台 Google Earth 
Engine (GEE) M MODIS 每 日 降雪 产品 MOD10A1 中 
分 别 逐 像 元 为 单位 识别 得 出 ,GEE 是 用 于 地 球 观测 
数据 分 析 的 云 平台 ”"。 用 户 可 以 在 https://earthen- 
gine.google.com 上 注册 后 访问 平台 。 在 该 平台 上 ， 
MODIOAL 产品 与 其 他 卫星 图 像 和 产品 (例如 MO- 
DIS „Landsat ,Sentinel—-1 和 Sentinel-2) 在 全 球 范 围 内 
被 使 用 并 经 常 更 新 。 如 果 栅 格 和 矢量 数据 不 直接 
可 用 ,也 可 以 很 容易 地 将 其 在 平台 中 提取 。 在 本 研 
究 中 ,使 用 GEE 编 程 语言 基于 Web 的 集成 开发 环境 
IDE 计算 积 雪 终日 参数 。 代 码 开 发 包括 数据 提取 、 
预 处 理 . 积 雪 终日 的 计算 统计 分 析 以 及 最 终 地 图 
的 导出 都 在 该 平台 上 进行 。 由 于 长 时 间 序 列 (新 疆 
日 数据 ) 的 计算 可 能 非常 密集 ,因此 采用 分 层 解 
决 方案 。 首 先 ,动态 数据 处 理 被 用 于 代码 开发 和 调 


不 能 超过 124 d 的 像素 。 掩 膜 冯 感 影像 上 的 水 , 仅 
分 析 在 陆 面 上 的 像素 。 最 终 经 过 以 上 3 种 处 理 后 再 
提取 和 统计 当年 1 月 1 日 开始 0% 被 雪 履 盖 的 最 后 
一 天 。 并 基于 行政 边界 剪裁 得 到 2000 一 2019 年 的 
新 疆 积 雪 终 日 数据 ,最 后 将 影像 导出 到 本 地 。 
(2) 一 元 线性 趋势 分 析 。 一 元 线性 回归 模型 可 
以 模拟 空间 上 每 个 像 元 的 积 雪 终日 时 间 变 化 ,斜率 
(Slope) 表 示 积 雪 终 日 时 间 变 化 的 趋势 ,用 最 小 二 乘 
法 来 计算 ,公式 如 下 : 
ndiixp)- Didp. 


Slope = 一 :一 
mn = 
式 中 :n 为 累计 年 数 ,n=20;p; 为 第 i 年 融雪 时 间 值 ; 
Slope 为 一 元 线性 回归 方程 的 斜率 ,是 n 年 连续 年 份 
的 趋势 , 当 Slope>0 表示 随时 间 i 的 增加 , 积 雪 终日 
处 于 推迟 趋势 ,反之 则 是 提前 趋势 。 

(3) 变异 系数 法 。 利 用 变异 系数 (Coefficient of 
variation,C,) 对 研究 区 2000—2019 年 积 雪 终 日 的 空 
间 格 局 变化 进行 稳定 性 分 析 ,C. 作 为 反映 观测 值 变 
异 程度 的 一 个 重要 指标 ,为 标准 差 与 多 年 平均 值 的 
比值 。 计 算 公 式 如 下 : 

c,= 2m (2) 
SM 
式 中 : Soy 为 每 个 栅 格 单元 的 标准 差 ; SM 为 积 雪 终 
日 的 多 年 平均 值 。 


3 结果 与 分 析 


3.1 新 疆 积 雪 终 日 的 空间 变化 特征 

根据 积 雪 终 日 的 分 布 情况 ,将 新 疆 分 为 北 疆 、 
南 疆 和 天 山 3 个 区 域 进行 积 雪 终 日 的 空间 分 析 。 如 
图 2 所 示 ,2000 一 2019 年 新 疆 积 雪 终 日 空间 分 布 存 


(1) 


试 。 其 次 ,将 每 年 和 一 年 中 不 同时 期 的 积 雪 产品 另 
存 ,然后 对 这 些 文件 进行 最 终 的 统计 分 析 。 

(1) 遥感 数据 处 理 。 基 于 GEE OE A, ,根据 影 
像 的 获取 时 间 影像 质量 ,快速 完成 时 间 序 列 遥 感 
数据 的 获取 和 预 处 理 ,包括 云 .阴影 掩 膜 .拼接 R 
剪 等。 本文 以 像素 为 单位 计算 每 年 的 积 雪 终日 ,最 
终 在 经 过 以 下 几 项 处 理 后 进行 统计 分 析 : 撼 膜 很 少 
被 雪 有 履 盖 的 像素 ,只 分 析 在 当年 1 月 1 日 到 12 月 31 
日 至 少 有 14d 被 10% 的 雪 和 覆盖 的 像素 。 掩 腊 长 时 
间 被 雪 有 覆盖 的 像素 ,在 一 年 内 10% 被 雪 履 盖 的 日 数 


在 明显 的 差异 ,整体 上 呈现 出 北 又 的 积 雪 终日 晚 于 
南 疆 ,山区 晚 于 分 地 , 且 两 者 的 空间 分 布 状况 存在 
明显 的 差异 。 北 弓 阿 尔 泰山 脉 一 带 、 仇 犁 河谷 以 及 
天 山西 段 .昆仑 山 一 带 积 雪 终 日 比 其 他 区 域 显 著 
晚 ,在 299~365 d 之 间 。 北 性 准噶尔 盆地 和 伊犁 河 
谷 积 雪 终日 明显 晚 于 南 疆 ,在 76~152 4 之 间 。 哈 密 
以 及 南 疆 塔里木 命 地 积 雪 终日 明显 早 于 北 疆 ,在 0~ 
61 d 之 间 。 

根据 积 雪 终 日 的 空间 分 布 和 变化 情况 可 以 看 
出 (图 2) , 随 着 气候 变 暖 北 狂 和 南 疆 .山区 和 盆地 积 
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(b) 20004 


(a) 年 平均 值 


积 雪 终日 /d 
61 122 183 244 305 365 


图 2 新 疆 积 雪 终 日 多 年 平均 与 年 际 变化 空间 分 布 


Fig.2 Spatial patterns of annual average and annual change in snow cover cessation day of Xinjiang 


雪 终 日 出 现 不 同 的 变化 。 南 疆 塔里木 分 地 和 东 吐 
区 域 积 雪 终日 从 2005 年 开始 逐渐 提前 , 北 疆 准噶尔 
贫 地 和 伊犁 河谷 逐渐 推迟 , 北 绪 和 南 疆 高 海拔 山区 
职 雪 终日 也 逐渐 推迟 ,由 此 可 见 在 不 同 的 海拔 高 度 
和 不 同 的 经 纬度 积 雪 终日 的 变化 有 所 不 同 。 
3.2 新 疆 积 雪 终 日 研究 期 初 终 对 比 变化 特征 

根据 积 雪 终日 数据 中 选择 2019 年 与 2001 年 2a 
数据 并 计算 出 两 者 之 差 , 比较 积 雪 终 日 初 终 的 差异 
(图 3)。 由 图 3 可 以 看 出 ,2019 年 和 2001 年 存在 明 
显 的 差异 ,整体 上 呈现 出 2019 年 南 性 塔里木 盆地 、 
阿尔 泰山 脉 一 带 以 及 东 疆 积 雪 终日 比 2001 年 提前 ， 
而 北 杜 准噶尔 盆地 天山 中段 和 昆仑 山脉 一 带 比 
2001 年 有 明显 的 推迟 。 北 又 准噶尔 盆地 和 南 疆 塔 
里 木 盆地 空间 分 布 状 况 恰 好 相反 ,塔里木 盆地 2019 
年 积 雪 终日 比 2001 年 提前 ,准噶尔 盆地 比 2001 年 推 
迟 。 北 疆 和 南 肚 ,高 海拔 山区 和 念 地 变化 特征 恰好 
相反 。 

2019 年 和 2001 年 新 疆 积 雪 终 日 差异 中 , 积 雪 终 
日 提前 的 区 域 占 新 疆 总 面积 的 26% , 积 雪 终日 推迟 
的 区 域 占 新 疆 总 面积 的 15%, 基 本 不 变 的 区 域 占 总 
面积 的 59%。 新 疆 积 雪 终日 提前 的 面积 大 于 推迟 


注 :红色 表示 积 雪 终 日 提前 , 蓝 色 表示 积 雪 终日 推迟 。 下 同 。 
图 3 新 疆 积 雪 终日 初 终 差 异 空间 分 布 


Fig. 3 Spatial pattern of the difference in snow cover 


cessation day of Xinjiang 


的 面积 。 
3.3 新 疆 积 雪 终 日 趋势 分 析 

图 4 表示 近 20 a 来 新 疆 积 雪 终日 变化 趋势 空间 
分 布 情 况 。 从 图 4 中 可 以 看 出 , 积 雪 终日 在 准噶尔 
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图 4 近 20 a 新 疆 积 


盆地 西部 和 南部 .阿尔 泰山 脉 一 带 、 天 山 山 脉 和 昆 
仑 山脉 部 分 高 海拔 区 域 有 推迟 趋势 ,最 大 推迟 幅度 
达到 16 d ,面积 占 新 疆 总 面积 的 8%。 准 踢 尔 盆地 中 
部 .伊犁 河谷 西部 和 塔里木 盆地 东部 基本 不 变 , 面 
积 占 新 疆 总 面积 的 48% ; 积 雪 终 日 在 塔里木 盆地 、 
东 疆 地 区 .准噶尔 盆地 北部 区 域 . 天 山南 坡 和 部 分 
高 海拔 区 域 有 提前 趋势 ,轻微 提前 和 明显 提前 区 域 
面积 约 占 新 疆 总 面积 的 44%。 积 雪 终 日 明显 提前 
地 区 主要 分 布 在 天 山南 坡 .塔里木 岔 地 西北 地 区 和 
东 疆 地 区 ,提前 幅度 达到 14 d。 积 雪 终 日 在 南 疆 塔 
里 木 命 地 和 东 疆 等 区 域 总 体 上 有 提前 趋势 , 北 疆 准 
噶 尔 盆地 积 雪 终日 呈现 提前 趋势 或 保持 不 变 。 在 
阿尔 泰山 脉 一 带 .天 出 山脉 中 部 .昆仑 山脉 等 区 域 ， 
积 雪 终日 总 体 上 呈现 提前 趋势 ,只 有 部 分 区 域 存在 
轻微 推迟 趋势 。 新 疆 积 雪 终 日 基本 不 变 的 区 域 面 
积 占 新 疆 总 面积 的 48% ,提前 趋势 的 区 域 面 积 占 新 
疆 总 面积 的 44% ,推迟 趋势 的 区 域 面积 占 新 疆 总 面 
积 的 8%。 新 疆 积 雪 终 日 呈现 提前 趋势 的 区 域 面积 
明显 大 于 推迟 趋势 的 区 域 面积 。 
3.4 积 雪 终日 年 际 波动 程度 

图 5 表示 经 波动 值 栅 格 计算 后 得 到 的 新 疆 近 
20 a 的 积 雪 终日 年 际 波动 程度 。 本 文 将 C. 值 定义 为 
显著 波动 ,轻微 波动 和 中 度 波动 3 个 波动 区 间 , 具 体 
分 类 情况 见 表 1。 

新 疆 积 雪 终日 整体 呈现 中 度 波动 居多 ,中 度 波 


雪 终 日 变化 趋势 


Fig.4 Variation trend of the snow cover cessation day in Xinjiang in recent 20 years 


动 和 显著 波动 接近 的 特征 。 积 雪 终 日 呈现 显著 波 
动 的 面积 占 新 疆 总 面积 的 26.12% ,主要 分 布 在 南 吐 
塔里木 盆地 中 部 和 南部 ;呈现 中 度 波 动 的 面积 占 总 
面积 的 40.31% ,主要 分 布 南 疆 ; 呈 现 轻微 波动 的 面 
积 占 总 面积 的 33.57% ,主要 分 布 在 北 疆 。 由 此 可 
见 ,新 疆 积 雪 终 日 近 20 a 来 不 稳定 , 南 疆 波动 性 
最 大 。 
3.5 ”新疆 积 雪 终日 的 年 平均 时 间 变 化 特征 
根据 积 雪 终 日 的 分 布 情况 ,将 新 疆 分 为 北 疆 、 
南 疆 和 天 山中 段 3 个 区 域 进行 积 雪 终 日 的 时 间 序 列 
分 析 ; 根 据 天 山中 段 . 北 疆 . 南 疆 积 雪 终日 数据 ,分 
别 计算 了 天 山中 有 段 , 北 疆 、 南 疆 区 域内 积 雪 终日 的 
年 平均 值 。 从 图 6 中 可 以 看 出 ,天 山中 段 区 域内 积 
雪 终 日 年 平均 值 2001 年 和 2004 年 开始 增加 ,2003 
年 和 2006 年 降低 ,天 山中 段 积 雪 终 日 随 年 际 变化 显 
著 , 呈 现 少 见 的 “M” 特 征 , 即 多 年 积 雪 年 份 中 明显 的 
波峰 和 波 谷 各 有 2 个 (图 6a) ,波峰 位 置 出 现在 2003 
年 和 2006 年 , 波 谷 位 置 则 出 现在 2004 年 和 2009 年 ， 
之 后 天 山中 段 积 雪 年 均值 呈现 不 稳定 状况 。 北 疆 
积 雪 终日 年 均值 从 2000 一 2006 年 基本 稳定 70~80 d 
之 间 ,2006 年 开始 出 现 不 稳定 状况 ,2006 一 2008 年 
积 雪 终日 年 均值 迅速 减少 ,之 后 2008 一 2010 年 迅速 
上 升 , 积 雪 终 日 年 均值 增加 , ASE 2008—2019 年 呈 
现 不 稳定 状况 (图 6b)。 根 据 南 疆 积 雪 终 日 数据 , 计 
算 了 南 疆 区 域内 积 雪 终 日 的 年 平均 值 ,从 结果 可 以 
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图 5 新 疆 年 际 积 雪 终 日 的 波动 程度 


Fig. 5 Variation coefficient of interannual snow cover cessation day in Xinjiang 


表 1 新 疆 年 际 积 雪 终日 的 C, 空 间 分 布 分 类 


Tab. 1 Classification of C, spatial distribution of interannual snow cover cessation day in Xinjiang 


类 型 波动 程度 C, 取 值 范围 像 元 数 面积 比 列 /% 
轻微 波动 0.037<C.< 1.038 2909778 33.57 
积 雪 终日 中 度 波动 1.038<C.< 2.320 3493170 40.31 
显著 波动 2.320<C< 4.514 2263823 26.12 


注 :C. 为 变异 系数 。 
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6 2000 一 2019 年 新 疆 积 雪 终日 年 平均 时 间 序 列 变化 


Fig. 6 Variations of the annual average snow cover cessation day in Xinjiang from 2000 to 2019 
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注 :红色 表示 集中 分 布 的 区 域 , 红 至 蓝 表示 集中 分 布 逐 渐 减 少 。 
7 2000 一 2019 年 新 疆 积 雪 终 


日 多 年 均值 与 海拔 高 度 的 关系 


Fig.7 Relationship between the multi-year mean value of snow cover cessation day and altitude in Xinjiang from 2000 to 2019 


看 出 , FASE 2003 年 开始 积 雪 终日 年 均值 增加 ,2006 
年 达到 最 高 值 之 后 逐渐 减少 2007 年 年 均值 最 小 ， 
2003 一 2009 年 呈现 不 稳定 状况 (图 6c)。 
3.6 新 疆 积 雪 终日 与 海拔 高 度 之 间 的 关系 

利用 DEM 数据 和 积 雪 终日 多 年 平均 数据 进项 
了 分 析 和 人 研究, 新疆 年 均 积 雪 终 日 与 海拔 高 度 之 间 
的 关系 中 (图 7) ,新疆 南 疆 、 北 疆 和 天 山中 段 积 雪 终 
日 在 不 同 的 海拔 高 度 上 具有 较 大 的 差异 性 。 北 性 
大 部 分 积 雪 终日 主要 分 布 在 海拔 高 度 为 130~1000 
m 的 区 域 , 积 雪 终日 为 60~90 d 之 间 , 积 雪 终日 随 着 
海拔 高 度 的 升 高 且 呈 增加 的 趋势 。 南 疆 大 部 分 积 
雪 终 日 主要 分 布 在 海拔 高 度 为 9000~1000 m 的 区 域 ， 
积 雪 终日 为 0~25 d 之 间 , 海 拔高 度 大 于 3000 m 的 一 


部 分 区 域 随 着 海拔 高 度 增加 积 雪 终 日 开始 增加 ,但 
大 部 分 高 海拔 区 域 在 0~20 d 之 间 , 主 要 跟 降水 气流 
来 源 密切 相关 。 天 山中 段 降 水 量 少 的 区 域 积 雪 终 
日 主要 分 布 在 海拔 高 度 为 500~1500 m 的 区 域 , 积 雪 
终日 是 0~30 4 之 间 。 天 山中 段 降水 量 高 区 随 着 海 
拔高 度 的 增加 积 雪 终日 天 数 增加 明显 , 积 雪 终 日 主 
要 出 现在 海拔 高 度 200~1500 m 区域 , 积 雪 终日 为 
75~100 d。 北 疆 积 雪 终 日 与 海拔 高 度 线 性 相关 性 较 
强 0.35, 其 次 是 天 山中 段 为 0.20, 最 低 的 南 疆 为 
0.18, 


4 讨论 


本 文 仅 结合 遥感 时 间 序 列 数 据 对 新 疆 积 雪 终 
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日 时 空 分 布 特征 进行 了 分 析 与 研究 ,今后 的 工作 中 
还 需要 着 重 将 积 雪 遥 感 数据 与 其 他 辅助 资料 相 结 
合 对 新 疆 积 雪 期 , 积 雪 厚度 的 变化 和 趋势 进行 深入 
人 研究。 此 外 , 积 雪 终日 年 均值 呈现 出 “M" 型 特征 的 
原因 及 其 形成 机 理 不 其 明度, 有待 今 后 进一步 深入 
人 研究。 本 研究 只 对 全 年 的 积 雪 终日 时 空 变化 特征 
进行 了 分 析 , 未 来 有 必要 增加 年 、 季 尺度 气象 内 子 
与 积 雪 终日 的 相关 性 研究。 同时 ,也 需要 采用 更 长 
时 间 序 列 的 积 雪 数据 ,对 新 疆 积 雪 终日 时 空 变化 规 
律 及 积 雪 时 空 变化 特征 与 气候 因子 之 间 的 响应 关 
系 进行 更 加 深入 的 研究 。 

新 疆 积 雪 终 日 时 空 变化 特征 中 南 疆 2004 一 
2005 年 呈现 不 稳定 状况 ,之 后 南 弓 积 雪 终 日 变化 较 
稳定 。 新 疆 积 雪 终 日 整体 上 呈现 出 的 不 稳定 状况 
可 能 与 新 疆 的 特殊 地 势 .特殊 地 貌 .人 研究 数据 采集 
时 间 序 列 有 关 。 在 某 种 程度 上 ,这 不 仅 反 映 了 新 疆 
积 雪 终日 年 内 变化 的 不 稳定 性 ,也 表现 出 不 同 区 域 
积 雪 变化 的 地 域 特殊 性 。 积 雪 终 日 年 均值 和 海拔 
高 度 之 间 的 关系 而 言 ,南星 、 北 疆 和 天 山中 段 积 雪 


雪 终 日 推迟 的 区 域 面积 占 新 疆 总 面积 15%, 主 要 分 
布 在 北 杜 准噶尔 盆地 、 天 山中 段 和 昆仑 山脉 一 带 。 
北 疆 准噶尔 盆地 天 山中 段 和 昆仑 山脉 一 带 积 雪 终 
日 2019 年 比 2001 年 有 明显 的 推迟 特征 。 

(3) 根据 2000 一 2019 年 新 疆 积 雪 终日 趋势 变 
化 来 看 , 积 雪 终日 提前 的 区 域 面积 占 新 疆 总 面积 的 
44% ,主要 分 布 在 南 疆 塔里木 合 地 和 东 疆 。 积 雪 终 
日 延迟 的 区 域 面积 占 新 疆 总 面积 的 8% ,主要 分 布 
在 阿尔 泰山 脉 .准噶尔 盆地 西部 .天 山中 段 . 北 对 和 
南 疆 高 海拔 山区 。 积 雪 终 日 不 变 的 区 域 面积 占 新 
性 总 面积 的 48% , 主要 分 布 在 准噶尔 盆地 中 部 。 南 
性 塔里木 岔 地 和 北 疆 准噶尔 盆地 有 明显 的 相反 趋 
势 。 从 年 际 波动 程度 上 看 ,新疆 积 雪 终 日 近 20 a 来 
旦 现 不 稳定 状况 ; 按 不 同 地 理 位 置 , 北 疆 轻 微波 动 ， 
南 疆 波动 性 最 大 , 东 疆 呈现 中 度 波 动 。 

(4) 时 间 上 ,2000 一 2019 年 北 疆 和 天 山中 段 积 
雪 终 日 年 均值 变化 较为 明显 ,天 山中 段 呈 现 “M” 型 
的 波动 特点 ,年 均值 在 2003 年 和 2006 年 处 于 高 位 
水 平 ,在 2009 年 处 于 低谷 水 平 。 北 疆 2006 年 开始 


终日 在 不 同 的 海拔 高 度 上 具有 较 大 的 差异 性 ,主要 
跟 降水 气流 来 源 密切 相关 。 在 全 球 气候 变 瞬 的 萌 
景 下 ,新 疆 气候 也 在 经 历 着 显著 变 上 暧 ,可 以 看 出 全 
球 气候 变 暧 对 新 疆 不 同 区 域 积 雪 终 日 的 影响 有 所 
不 同 。 


5 结论 


(1) 新 疆 2000—2019 年 积 雪 终日 具有 明显 的 
空间 分 异 特征 , 积 雪 终日 空间 分 布 体现 北 疆 大 于 南 
吐 、 东 部 大 于 西部 的 分 布 特征 。 北 疆 阿 尔 泰山 一 
带 .伊犁 河谷 和 准噶尔 盆地 以 及 天 山中 段 和 南 性 高 
海拔 区 域 积 雪 终 日 较 长 ,哈密 以 及 南 性 塔里木 盆地 
低 海拔 区 域 积 雪 终日 较 短 , 北 疆 的 积 雪 终日 分 布 明 
显 长 于 南 疆 的 积 雪 终 日 。 根 据 积 雪 终日 多 年 变化 
可 以 看 出 ,从 2005 年 开始 南 疆 塔里木 倪 地 积 雪 终日 
逐渐 缩短 , 北 疆 积 雪 终 日 逐渐 变 长 ,空间 分 布 与 多 
年 情况 有 明显 的 变化 。 

(2) 根据 2019 年 和 2001 年 新 疆 积 雪 终 日 年 度 
差异 变化 特征 中 可 以 看 出 :2019 年 和 2001 年 积 雪 终 
日 提前 的 区 域 面 积 占 新 疆 总 面积 的 26% , 主要 分 布 
在 南 疆 塔里木 盆地 和 阿尔 泰山 脉 一 带 和 东 疆 区 
域 。 基 本 不 变 的 区 域 面积 占 新 疆 总 面积 的 59% , 积 


呈现 不 稳定 状况 ,2008 一 2015 年 积 雪 终 日 年 均值 迅 
速 减 少 处 于 低谷 水 平 ,2010 年 积 雪 终 日 年 均值 处 于 
高 位 水 平 。 北 性 随 着 海拔 高 度 增 加 积 雪 终日 变 长 ， 
南 疆 和 天 山中 段 积 雪 终 日 与 海拔 高 度 之 间 的 相关 
性 较 复杂 。 


参考 文献 (References) 


[1] Barnhart T B, Molotch N P, Livneh B, et al. Snowmelt rate dictates 
streamflow[J]. Geophysical Research Letters, 2016, 43(15): 8006- 
8016. 

[2] Bengtsson L. Snowmelt estimated from energy budget studies|J]. 
Hydrology Research, 1976, 7(1): 3-18. 

[3] Duchesne L, Houle D, D’ Orangeville L. Influence of climate on 
seasonal patterns of stem increment of balsam fir in a boreal forest 
of Quebec, Canada[J]. Agricultural and Forest Meteorology, 2012, 
162: 108-114. 

[4] Kudo G, Nordenhäll U, Molau U. Effects of snowmelt timing on 
leaf traits, leaf production, and shoot growth of alpine plants: Com- 
parisons along a snowmelt gradient in northern Sweden[J]. Coence, 
1999, 6(3): 439-450. 

[5] Solomon S, Qin D, Manning M, et al. Climate change 2007: The 
physical science basis. Contribution of working group I to the 
fourth assessment report of the intergovernmental panel on climate 
change (IPCC)[J]. Computational Geometry, 2007, 18(2): 95-123. 


[6] Beniston M. Impacts of climate change on water and associated 


202208.00061v1 


chinaXiv 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


买 合 木 提 江 : 维 吉 日 等 : 


economic activities in the Swiss Alps[J]. Journal of Hydrology, 
2010, 412: 291-296. 

Clow D W. Changes in the timing of snowmelt and streamflow in 
Colorado: A response to recent warming[J]. Journal of Climate, 
2010, 23(9): 2293-2306. 

Harpold A A, Brooks P D, Rajogopalan S, et al. Changes in snow- 
pack volume and snowmelt timing in the Intermountain West[J]. 
Water Resources Research, 2012, 48(11): 1-11. 

魏 月 . JEJE K J Es REZE FER TNE FED]. 马 鲁 木 齐 : 
新 疆 师 范 大 学 , 2010. [Wei Yue. Remote sensing monitoring of 


snow depth change in north part of Xinjiang[D]. Urumqi: Xinjiang 
Normal University, 2010. ] 

TAH ohare. 天 山 积 雪 初步 研究 由 . 干旱 区 地 理 , 1992, 15 
(3): 9-21. [Qiu Jiaqi, Sun Xihua. Preliminary study on snow cover 
in the Tianshan Mountain[J]. Arid Land Geography, 1992, 15(3): 
9-21. ] 

张 一 驰 , 李 宝 林 , 包 安 明 , 等 . IARI fat a aS A EL]. 
中 国 科 学 : 地 球 科学 , 2006, 36( 增 刊 2): 24-32. [Zhang Yichi, Li 
Baolin, Bao Anming, et al. Simulation study on snowmelt runoff in 
Kaidu River Basin[J]. Scientia Sinica (Terrae), 2006, 36(Suppl. 2): 
24-32. | 

杨 绍 富 . Fi LEP KP K FEV), ESE Ti 8) 
[D]. 乌鲁木齐 : 新 疆 大 学 , 2009. [Yang Shaofu. The study of wa- 


ter and heat coupling experimentin snow melt processing: Take 
Juntang River Basin Watershed in northern slope of Tianshan 
Mountains as an example[D]. Urumqi: Xinjiang University, 2009. | 
岳 春 芳 , 何 训 江 , SCA TH BM «SEB, 金 沟 河流 域 水 权 分 配 及 转 
让 研究 加 . 水 资源 与 水 工程 学 报 , 2012, 23(3): 48-50. [Yue 
Chunfang, He Xunjiang, Aisha Ailimiguli. Research on the alloca- 


tion and transfer of water rights in Jingouhe River Basin[J]. Journal 
of Water Resources and Water Engineering, 2012, 23(3): 48-50. ] 
Qin Y, Abatzoglou J T, Siebert S, et al. Agricultural risks from 
changing snowmelt|J]. Nature Climate Change, 2020, 10(5): 459- 
465. 

Diffenbaugh N S, Scherer M, Ashfaq M. Response of snow-depen- 
dent hydrologic extremes to continued global warming[J]. Nature 
Climate Change, 2013, 3(4): 379-384. 

唐 尔 泗 , 程 国 栋 , 董 增 川 . 中 国 西北 干旱 区 冰雪 水 资源 与 出 山 径 
流 [M]. 北京 : 科学 出 版 社 , 2002: 1-2. [Tang Ersi, Cheng Guodong, 


lly 


Dong Zengchuan. Snow and ice water resources and mountain run- 
off in arid areas of northwest China[M]. Beijing: Science Press, 
2002: 1-2. ] 

王 秋 香 , 魏 文 寿 , TER. 新 疆 北 疆 最 大 积 雪 深度 EOF 展开 场 
的 时 间 变 化 规律 [. 冰川 冻 土 , 2008, 30(2): 244-249. [Wang 


Qiuxiang, Wei Wenshou, Wang Jinmin. Temporal variation of the 


maximum snow cover depth in north Xinjiang derive from EOF\J]. 
Journal of Glaciology and Geocryoiogy, 2008, 30(2): 244-249. | 
Guo L P, Li L H. Variation of the proportion of precipitation occur- 


ring as snow in the Tianshan Mountains, ChinalJ]. International 


2000—2019 4] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


FE 新 疆 积 雪 终日 时 空 变化 特征 


ChinaXiv 合 作 期 刊 


T 


Journal of Climatology, 2015, 35(7): 1379-1393. 

ett, 丁 建 丽 , BORA, 等 . 2001 一 2015 年 天 山 山区 积 雪 时 空 变 
化 及 其 与 温度 和 降水 的 关系 町 . 冰川 冻 土 , 2018, 40(2): 249- 
260. [Qin Yan, Ding Jianli, Zhao Qiudong, et al. Spatial-temporal 


variation of snow cover in the Tianshan Mountains from 2001 to 
2015, and its relation to temperature and precipitation[J]. Journal 
of Glaciology and Geocryology, 2018, 40(2): 249-260. | 

王国 亚 , GARIR, SOR, 等 . 新 疆 阿勒泰 地 区 积 雪 变化 特征 及 其 
对 冻 土 的 影响 四. 冰川 冻 土 , 2012, 346): 1293-1300. [Wang Guoya, 


Mao Weiyi, He Bin, et al. Changes in snow covers during 1961— 
2010 and its effects on frozen ground in Altay region, Xinjiang[J]. 
Journal of Glaciology and Geocryology, 2012, 34(6): 1293-1300. ] 

DA, 刘 时 银 , REE, 等 . 2002—2018 年 叶 尔 羌 河流 域 积 雪 时 空 
变化 研究 [J]. 干旱 区 地 理 , 2021, 44(1): 15-26. [Yi Ying, Liu Shi 


yin, Zhu Yu, et al. Spatiotemporal variation of snow cover in the 
Yarkant River Basin during 2002—2018[J]. Arid Land Geogra- 
phy, 2021, 44(1): 15-26. ] 

UE Ma, 李 忠 勤 , KE. 2001—2017 年 祁连山 积 雪 面 积 时 空 变 
化 特征 吓 . 干旱 区 地 理 , 2019, 42(1): 56-66. [Liang Pengbin, Li 


Zhongqin, Zhang Hui. Temporal-spatial variation characteristics of 
snow cover in Qilian Mountains from 2001 to 2017[J]. Arid Land 
Geography, 2019, 42(1): 56-66. ] 

WACK, 苏 宏 超 , 王国 亚 , 等 . DEDKI ARORA AC 
M: 水 文 效应 由. 冰川 冻 土 , 2013, 35(3): 513-527. [Shen Yong- 


ping, Su Hongchao, Wang Guoya, et al. The responses of glaciers 
and snow cover to climate change in Xinjiang (I): Hydrological ef- 
fect[J]. Journal of Glaciology and Geocryology, 2013, 35(3): 513- 
527.] 

刘 吉 强 . 基于 空间 分 析 法 的 新 疆 融 雪 洪 水 径流 模拟 叫 . 黑龙 江 
水 利 科 技 , 2020, 48(10): 17-20. [Liu Jiqiang. Runoff simulation 


of snowmelt flood in Xinjiang based on spatial analysis method[J]. 
Heilongjiang Water Conservancy Science and Technology, 2020, 
48(10): 17-20. ] 

杨 金 明 , 李 诚 志 , 房 世 峰 , 等 . 新 疆 地 区 季节 性 融雪 洪水 模拟 与 
THAR OF SE (J). 新 疆 大 学 学 报 (自然 科学 版 ), 2019, 36(1): 80-88. 
[Yang Jinming, Li Chengzhi, Fang Shifeng, et al. Simulation and 


forecast of seasonal snowmelt flood in Xinjiang[J]. Journal of Xinji- 
ang University (Natural Science Edition), 2019, 36(1): 80-88. | 

王 冉 旋 . 基于 气象 数值 预报 的 高 精度 融雪 径流 预报 技术 研究 
U]. 科技 创新 导报 , 2019, 16(32): 21-23. [Wang Ranxuan. Research 


on high-precision snowmelt runoff prediction technology based on me- 
teorological numerical prediction[J]. Science and Technology Inno- 
vation Herald, 2019, 16(32): 21-23. ] 

张 云 云 , 张 航 涛 , 师 庆 东 , 等 . 天 山北 坡 中 段 草 地 , 林地 积 雪 消融 
过 程 的 定量 化 分 析 町 . 水 土 保持 学 报 , 2019, 33(3): 108-114. 
[Zhang Yunyun, Zhang Yutao, Shi Qingdong, et al. Quantitative 


analysis of snow melting process in grassland and forest land on 
the middle section of northern slope of Tianshan Mountain[J]. Jour- 


nal of Soil and Water Conservation, 2019, 33(3): 108-114. | 


202208.00061v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


FHEWY 


[28] 梁 建 辉 .气候 变化 环境 下 新 疆 喀什 冰川 河流 融雪 径流 动态 响应 [30] AER, 黄 韶华 . 对 新 疆 生 态 环境 问题 及 其 对 策 的 若干 思考 四]. 
区 


研究 中. 水 利 规划 与 设计 , 2017(5): 41-43, 62. [Liang Jianhui. 干旱 区 资源 与 环境 , 1999(4): 1-8. [Zhou Huarong, Huang Shao- 
Study on dynamic response of snowmelt runoff of Xinjiang Kash- hua. On Xinjiang ecological enviromental problem and its counter- 
gar glacier river under climate change environment|J]. Water Re- measures[J]. Journal of Arid Land Resources and Environment, 
sources Planning and Design, 2017(5): 41-43, 62. ] 1999(4): 1-8. ] 

[29] 李 培 基 . 中 国 季节 积 雪 资源 的 初步 评价 由. 地 理学 报 , 1988(2): [31] Gorelick N, Hancher M, Dixon M, et al. Google Earth Engine: 
108-119. [Li Peiji. Preliminary evaluation of seasonal snow resourc- Planetary-scale geospatial analysis for everyone[J]. Remote Sens- 
es in China[J]. Acta Geographica Sinica, 1988(2): 108-119. ] ing of Environment, 2017, 202: 18-27. 


Temporal and spatial changes of the snow cover cessation day from 2000 to 
2019 in Xinjiang, China 


Maihemutijiang WEIJIDAN', Yusufujiang RUSULI'*, QIU Zhongli 
(1. College of Geographic Science and Tourism, Xinjiang Normal University, Urumqi 830054, Xinjiang, China; 
2. Xinjiang Key Laboratory of Lake Environment and Resources in Arid Zone, Xinjiang Normal University, Urumqi 830054, 
Xinjiang, China) 


Abstract: In this study, based on Google Earth Engine (GEE) remote sensing cloud platform, using MODIS 
snow product from 2000 to 2019, the diurnal snow cover cessation day information in Xinjiang, China is 
extracted and calculated, and the diurnal spatial and temporal variation characteristics and variation trend of snow 
cover cessation day in Xinjiang are analyzed by using trend analysis, coefficient of variation methods. The results 
show the following: (1) with the Tianshan Mountains as the dividing line, and the snow cover cessation day of the 
north of the Tianshan Mountains is longer than that of the south. The mountainous area is the high-value area of 
snow cover cessation day, and the basin is the low value area of snow cover cessation day. The snow cover 
cessation day in the Junggar Basin and Ili River Valley in northern Xinjiang lasts 75-114 d all year long. The Tarim 
Basin in southern Xinjiang is a low-value area between 0 and 31 d. The Altai Mountains, Tianshan Mountains and 
Kunlun Mountains belong to high value areas between 224 and 365 days. (2) There are obvious spatio-temporal 
differences in snow cover cessation days in southern Xinjiang and northern Xinjiang. From 2000 to 2019, the 
snow cover cessation day in the Junggar Basin and the high-altitude mountains in northern Xinjiang has an 
obvious trend of delaying, with a delay of 14 days, accounting for 8% of the total area of Xinjiang. The Tarim 
Basin and the areas across eastern Xinjiang showed an obvious advance advancing trend, the advance advancing 
rate reached 16 d, accounting for 44% of the total area of Xinjiang. Tarim Basin and Junggar Basin have opposite 
trend. (3) The snow cover cessation day of Xinjiang varied all year, and central Tianshan Mountains and northern 
Xinjiang showed instability in changes of snow cover cessation day all year long. The snow cover cessation day 
of central Tianshan Mountains showed overall “M” characteristics shape curve, with the annual average snow 
melting days showing obvious peaks and troughs during 2002—2009. Moreover, northern Xinjiang showed large 
interannual changes with instability during 2009—2019 and with obvious peaks and troughs as well as large 
interannual variability. 


Key words: GEE cloud computing; snow cover cessation day; temporal and spatial distribution; change trend 


